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Abstract. It is shown theoretically that, in the superconducting state with d x 2_ y 2 -pairing, a 
strong Pauli paramagnetic depairing (PD) induces not only the modulated Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) superconducting state but also an incommensurate antiferromagnetic 
(AFM) (or spin-density wave) ordering with Q-vector nearly parallel to a gap node. In this 
mechanism of field-induced coexistence of the d-wave superconducting and AFM orders, a pair- 
density wave does not have to be assumed. It is argued that this is the common origin of both 
the coexistent FFLO and AFM phases of CeCoIns and the AFM quantum critical behavior 
around the superconducting H C 2{0) seen in several unconventional superconductors. 



1. Introduction 

At present, we have two intriguing issues on a magnetic order or fluctuation occurring close 
to H C 2(0) in (i-wave paired superconducting states. One is the antiferromagnetic (AFM) 
quantum critical behavior reflected in transport measurements around H C 2(0) of heavy fermion 
superconductors CeCoIns [1-3], pressured CeRhlns [4], NpPdsA^ [5], and Tl-based cuprates 
[6]. The other is the AFM order [7] in the high field and low temperature (HFLT) phase [8] of 
CeCoIns which has been identified, based on measurements [9,10] and theoretical explanations 
[11,12] of elastic properties and doping effects, with a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state [8,13]. It is notable that the materials listed above have a couple of common features, such 
as the ci-wave pairing and a strong Pauli paramagnetic depairing (PD). 

In this report, we point out that, in nodal d-wave superconductors, a field-induced 
enhancement of PD tends to induce an AFM ordering in the superconducting phase just below 
H C 2(0). Conventionally, the AFM ordering is suppressed by the superconducting ordering in 
zero field [14], and, even in nonzero fields, the quasiparticle damping effect brought by the AFM 
fluctuation usually suppresses a relatively weak PD effect, suggesting a competition between 
the two orderings. It is found, however, that a strong PD rather favors coexistence of a (i-wave 
superconductivity and an incommensurate AFM order, leading to an enhancement of AFM 
ordering or fluctuation just below H&ifi). 



2. Model and Calculation 

We work in a BCS-like electronic Hamiltonian [14] with the superconducting energy gap A and 
the AFM moment m. By treating A at the mean field level and m = |m| as a fluctuation, 



respectively, the free energy describing the two possible orderings in zero field is described by 
J-(H = 0) = |a! 3 r« ? - 1 |A(r)| 2 -rinTr CiC t im exp[-(^ Am - / uiV)/T], 

H Am -fiN = ^i|m(q)| 2 + ]T 4,a £( k ) <W c k ,/3 
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where *(q) = -i(a y ) a ,p Ek w k ^ k+q/2ia c k+q/2 ^/2, S„(q) = {p v ) a ,p Ek 4-q !( A+Q,/9/ 2 > 4,a 
creates a quasiparticle with spin index a and momentum k, a v are the Pauli matrices, \i is 
the chemical potential, and the positive parameters g and U are the attractive and repulsive 
interaction strengths leading to the superconducting and AFM orderings, respectively. The 
gap function % satisfies Wk+Q = — for the d x 2_ y 2 -pairing state, and the dispersion £(k) 
satisfies £(k) = — £(k + Q) + T c 5ic, where Q is the commensurate AFM modulation wavevector 
and is (tt, it) for the d x 2 _ y 2 -pairing. A small deviation from the perfect nesting is measured 
by a small parameter Sic f° r a nearly-free electron model. If the tight-binding model with 
dispersion £(k) = — ti(cos(k x a) + cos(k y a)) — t2Cos(k x a)cos(k y a) — \i and the lattice constant 
a is directly used in examining an AFM ordering, the corresponding incommensurability is 
measured by the second term of the above £(k). In a nonzero field (H ^ 0) of our interest, the 
Zeeman term jBH(a z ) a ^ needs to be added to £(k)<5 aj/ 3. At least in the case with a continuous 
-ff C 2-transition like Fig. 2 (a) below, the orbital depairing needs to be incorporated through the 
familiar quasiclassical treatment on the quasiparticle Green's function [13]. 

To see the position of the AFM ordering, it is convenient to examine the Gaussian AFM 
fluctuation term T m in the free energy J 7 , T m = TJ2n lndet[£/ _1 £q !q ' — Xq,q'(^) ]> where 

Xq , q ,(fi) = f T dr(T T SUcL;r)SM^))e iUT , (2) 
jo 

and S v (q;T) denotes S u (q) at imaginary time r. For the moment, we focus on the Pauli limit 
with no orbital depairing and with uniform A in which Xq,q' (^) = [x^ (<1> ^) + X^ an ^ ^)]^q,q' ; 
and T m = -TE qi! ilnI(q,fi), where X^faQ) = U' 1 '- X (n) (q,^) - x {an) (q,^)- A second 
order AFM ordering occurs when Xq 1 = X -1 (0, 0) = 0. The 0(|A| 2 ) terms in x an d X^ are 
expressed by Fig.l (a) and (b), respectively. They have been studied previously [14] in H = 
case, where x«(A) = x (n) (0,0) - X (n) (0, 0)|a=o + X (an) (0,0) taking the form 




(a) (b) 

Figure 1. Diagrams describing (a) x^ an d (b) X^ up to 0(|A| 2 ), where the cross denotes the 
particle-hole vertex on the AFM fluctuation, while the filled circle implies the particle-particle 
vertex on A or A*. 
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(3) 



behaves like T~ 2 in T —> limit and is negative so that the AFM ordering is suppressed by 
superconductivity [14]. In eq.(3), Ge,a(p) = (is - £p+7bH<j) _1 denotes the quasiparticle Green's 
function defined in the normal state. The last term implies x^ n \^i 0)- 

To explain effects of strong PD, let us first explain the m || H case in which a = a. In this 
case, the two terms in eq.(3) are found to take the same form as the coefficient of the 0(|A| 4 ) 
term of the superconducting Ginzburg-Landau (GL) free energy and thus, change their sign upon 
cooling [13]. Hence, Xs(A-) becomes positive for stronger PD, leading to a lower T m , i.e., an 
enhancement of the AFM ordering in the superconducting phase. As well as the corresponding 
PD-induced sign-change of the 0(|A| 4 ) term which leads to the first order //^-transition [13], 
the PD-induced positive Xs is also unaffected by inclusion of the orbital depairing. 

In m _L H where a = —a, a different type of PD-induced AFM ordering occurs in a d- 
wave pairing case with a gap node along Q where w p+ q = —w p : In this case, the first 
term of eq.(3) arising from x( n )(0,0) remains negative as in zero field case and becomes 
-iV(0)|A| 2 /[2(7B-ff) 2 ] in T ->• limit with no PD-induced sign change, where iV(0) is the 
normal density of states. Instead, the last term of eq.(3) implying x^ an -*(0, 0) and thus, Xs are 
divergent like N(0)[ \A\/(^bH) ] 2 |ln[Max(i, |<Sic|)]| in T -> limit while keeping their positive 
signs owing to the relation tw p +Q • w p < 0, where t = T/T c . This divergence is unaffected by 
including the orbital depairing. That is, in the d x 2_ y 2-w&ve case with Q = (ir, it), the AFM 
order tends to occur upon cooling in m _L H. In contrast, x^ an ^(0,0) is also negative in the 
d xy -w&ve case with the same Q satisfying w p w p+ q > so that the AFM ordering is suppressed 
with increasing H. 

3. Examples of Phase Diagrams 

In this section, examples of the resulting low temperature phase diagram near H C 2(0) will 
be presented. In the BCS-like model (1) and up to the 0(|A| 2 ) terms (see Fig.l), the U c i- 
transition, i.e., the mean field superconducting transition in H ^ 0, is of second order even 
at lower temperatures for a = 7si? or b(0)/(27rT c ) < 0.3 (see Fig.2(a)), where a is nothing but 
the Maki parameter except a difference in the numerical factor, while it becomes of first order 
for larger a ~ 1.1 [13]. It is reasonable to expect the former to correspond to the case of 
CeRhlns under a pressure [4]. Figure 2(a) is one of the phase diagrams in such a case, where 
the Neel temperature Tn in the normal state with perfect nesting or U was assumed to be the 
only parameter measuring the pressure dependence. The actual AFM transition temperature in 
H > H C 2(0) for Tn/T c = 0.02 and 0.35 are zero and less than 0.35T C , respectively, because of 
the finite <5ic — 0.6 used in the calculation. Reflecting the AFM ordering enhanced by PD, the 
decrease of Tn, corresponding to an increase of pressure, results in the shrinkage of the AFM 
phase just below H C 2(0), which reduces to an apparent AFM quantum critical point by a further 
increase of pressure. 

Figure 2(b) is the corresponding result in the tight binding model in the Pauli limit with 
no orbital depairing (vortices). The i/ C 2-transition is of first order in the temperature range 
shown there. Due to the discontinuous nature of the i/ C 2-transition, an apparent AFM quantum 
critical point is estimated, in h = H/H C 2(0) > 1, to lie at a lower field than H C 2(0) in spite of 
the PD-induced AFM ordering just below H C 2(0). This is consistent with the observations in 
CeCoIns [2,7]. We note that the anomalous doping effect in CeCoIns [10] cannot be explained 
without a spatial modulation of |A| in the HFLT phase [12], implying that both the AFM and 
FFLO orders coexist in the HFLT phase of CeCoIns. Calculation results in the case including 
the FFLO structure will be reported elsewhere [15]. 

We note that, in the present theory explaining the AFM order just below H C 2(0) in CeCoIns 
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Figure 2. Typical t (= T/T c ) v.s. h = H/H C 2(0) phase diagrams (a) following from the use of 
a = 0.3 leading to a second order iif C 2-transition even in low t limit and (b) in the Pauli limit 
with a first order /^-transition in low t, respectively. In both figures, an AFM phase can occur 
below a solid curve on which Xq 1 = 0, and each nearly vertical dotted curve is the corresponding 
if C 2(T)-curve. Note that, in Fig. 2(a), two AFM phase boundaries for T^/T c = 0.35 and 0.02 are 
shown in a single figure. Figure 2 (b) was obtained from the corresponding calculation in the 
tight binding model with the parameters U = 33, t\ = 100, and ti = 0.25 in the unit of T c and 
by taking account of the full A dependence with no limitation to the 0(|A| 2 ) term. The dashed 
curve denotes the possible upper limit of the AFM transition temperatures at which Xs = 0. 
The lower panel of (b) is the /i-dependence of X^ 1 at t = 0.05. 




with strong PD, the assumption [16] of an additional pairing channel (pair-density wave) is 
unnecessary, and that both of the AFM order [7] and other observations, such as the anomalous 
doping effect [10], in the HFLT phase of CeCoIns are explained consistently if the FFLO 
modulation in the HFLT phase is assumed. 

This work was financially supported by Grant-in-Aid for Scientific Research [No. 20102008 
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